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The B-B mixing, B-meson lifetimes, the Bs-Bs lifetime difference and SUSY FCNC effects in AB — 2 processes 
are very important measurable quantities in B-meson phenomenology whose theoretical predictions depend on un- 
known matrix elements of several four-fermion operators. We present preliminary results for the matrix elements 
of the relevant four-fermion operators computed on a sample of 600 lattices of size 24=* X 40 at /3 = 6.0, using the 
SW-Clover action for light quarks with rotated light quark propagators and the lattice version of the IfQET for 
heavy quarks. As a necessary ingredient of the calculation, we also present results for the next-to-leading order 
matching of the full theory to the lattice HQET (one- and two-loop anomalous dimensions, one-loop QCD-HQET 
matching coefficients and one-loop continuum-lattice HQET matching coefficients). 



1. Introduction and Motivation. 

B-hadron decays are a very important source 
of information on the physics of the standard 
model (SM) and beyond. In many important 
cases, long distance strong contributions to these 
processes can be separated into matrix elements 
of local operators. Lattice QCD can then be used 
to compute these non-perturbative parameters 
from first principles. A list of some four-fermion 
operators relevant to _B-meson phenomenology is 
the following: 



SUSY effects in AB = 2 transitions y, and 
QtLL and Qtfp; to spectator effects in the 
B-meson lifetimes 

Our aim is to non-perturbatively evaluate the 
matrix elements between S-meson states of all 
the operators relevant to phenomenology using 
lattice simulations of the HQET for the 6-quark. 
The procedure to perform the transition from 
QCD to lattice HQET, a combination of analytic 
and numerical calculations, consists in the four 
steps described in detail below. 



• the operator Oll = ^^7^ (1 - lb)qb^^ (1 - 
75)9 ( 9 denotes a light quark u, d or 
s) , which determines the theoretical predic- 
tion of the B°-B° mixing in the SM, to- 
gether with Ofi(^fl) = h{lT 15) qb{\T 
75) Q^ Orb. and Of ^j, with obvious notation, 
parametrize SUSY effects in Ai? = 2 tran- 
sitions [Q. Furthermore, Oll and Of^ de- 
termine the Bs width difference AF^, ||^; 

• Qll = 67^ (1-75)997^ (1-75) together 
with Qf ^, parametrize spectator effects in 
the tb and tb, lifetimes Q; 

• finally, the corresponding operators with a 
t"^ (the generator of the SU(3) group) inser- 
tion: OtLL, Otfj^ and Otf^ contribute to 



2. First step: the continuum QCD — 
HQET matching. 

The starting point is to express the QCD op- 
erators as linear combinations of the HQET ones 
in the continuum at a given scale, say, fj, = mi,. 
To do so, some amplitudes of the relevant opera- 
tors between appropriate external _B-meson states 
are evaluated to one-loop order both in QCD and 
HQET. Once the amplitudes have been renormal- 
ized using the same scheme in both theories, the 
QCD amplitud is expanded in powers of l/mi, 
to lowest order. Finally, the coefficients of the 
matching can be obtain by subtraction at the 
scale ^. It should be stressed that the matching 
is renormalization scheme dependent. Notice also 
that new effective operators can be generated. 

For Oll, the matching was determined by 
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Flynn et al using NDR ^ and by Gimenez ||^ 
using DRED. Now, we have extended the calcu- 
lation to all the relevant operators in section 1. 
As we said before, new effective operators arise at 
one-loop order; i.e. the operator Of is expressed 
as a linear combination of two HQET operators 
and Oll- We refer the reader to ref.[|j for 
details. 

3. Second step: running down to fi — a^^ 
in the HQET. 

The HQET operators obtained in step 1 at 
the high scale fi = uii, are evolved down to a 
lower scale /x = a~^, appropriate for lattice sim- 
ulations, using the HQET renormalization group 
equations. The running is governed by the cor- 
responding anomalous dimension matrices. At 
NLO, one has to calculate both the one- and two- 
loop anomalous dimension matrices of the rele- 
vant effective operators. The former are renor- 
malization scheme independent but the latter are 
not. The scheme dependence of the two-loop 
anomalous dimension matrices cancels with the 
scheme dependence of the matching in step 1. 

The calculation of the one-loop anomalous di- 
mension matrix is not difficult 
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in the basis {Oll, Of^, Olr, O^l)- 

However, only the two-loop anomalous dimen- 
sion of the operator Oll is known. It was de- 
termined by Gimenez § (see also ref.|0]) using 
NDR and DRED. Now we have extended the cal- 
culation to all effective operators in the basis. 
Technically this is a very hard task because it 
involves the evaluation of 49 two-loop Feynman 
diagrams in the HQET. Moreover, the operators 
mix among them, see eq.(|^), and with evanes- 
cent operators under renormalization. However, 
some useful relations between the matrix ele- 
ments of the anomalous dimension matrix can be 
obtained and used as a check of the calculation 
using a renormalization prescription preserving 
Fierz symmetry. We skip all details and refer the 



reader to ref.Q. 

4. Third step: continuum lattice HQET 
matching. 

Having obtained the continuum HQET opera- 
tors at the scale fi = a~^, they are expressed as 
a linear combination of lattice HQET operators 
at this scale. The procedure is very similar to 
step 1: one matches two amplitudes at one- loop 
order, one in the continuum HQET and the other 
in the lattice HQET. However, some subtle points 
to remember are the following: 

1. the matching has to be calculated using the 
same action used in the numerical simula- 
tion; the matching coefhcients depend on 
the action used in the simulation. 

2. Due to the breaking of chiral symmetry in- 
duced by the Wilson term for light quarks, 
the original operators can mix with new lat- 
tice operators. 

3. The matching, as in step 1, depends on 
the renormalization procedure used to de- 
fine the operators in the continuum HQET, 
which cancels with the remaining scheme 
dependence of the running in step 2. 

The lattice counterpart of the operator Oll is 
already known. It was determined by Flynn et al 
Q for the Wilson action and by Borrelli and Pit- 
tori H for the SW-Clover action. A gain, we have 
extended the calculation to all the HQET oper- 
ators in the basis of eq.(|l|). Our results for Oll 
agree with ref.Q but disagree with ref. [0 in the 
sign of the contribution proportional to O}^** (see 
below) of the group of Feynman diagrams with ra- 
diative corrections between light quark legs only. 
We refer the reader to ref. ^ for details. 

As we said before, new lattice operators arise 

Olatt r\latt , /^latt , n / y^latt S i r^latt S \ 

N — ^LR + ^RL + ^ I ^LR + ^ RL ) 

oT = 3/2(Oit* + O^^Ij*) 

+ 4(Oit^ + (2) 

Olatt r^latt , r\latt , n ( y^latt S i r^latt S \ 

P — ^LR. + <^i?,L + " I <^Li?, + ^RL ) 
Qlatt ^ Qlatt _^ QlM ^ g Qlatt S ^ Qlatt S -j 
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Figure 1. The ratios Rqll ? -^of ^ ? Ron ^'^^ ^Oa 
computed at K = 0.1432. 
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Figure 2. The ratios Rqp , Roq , Rolr '^^'^ -^of 
computed at K — 0.1432. 



We found that 0^"f mixes with 0^)j* and Oj^", 
O^''^* with O^f ^ and O^", O^"^*'^ with O^"^', 



and O^" and finally 0'^£f} mixes with O^** 

for de- 



and Og". We refer the reader to ref. | 
tails. 



5. Fourth step: lattice computation of the 
matrix elements. 

In order to obtain the B-parameters of the rel- 
evant operators, we simulate on the lattice the 
ratio (see, for example, ref.fiot]) 



Coi {—tiih) 
8/3C{-ti)C{t2) 



(3) 



where Co; (^1,^2) is the three-point correlation 
function with an insertion of the lattice oper- 
ator Oi and C{t) is the two-point correlation 
function for the S-meson. To extract the B- 
parameter, we search for a plateau in ti at fixed 
(and large) t2- We observe, see figs.l and 2, good 
plateaux over large time-distances for all oper- 
ators. This makes us confident that the light- 
est meson state has been isolated. Furthermore, 
our results are almost independent of the light 
quark mass, therefore they can be safely extrap- 
olated to the chiral limit. From figs.l and 2, 
we see that the correction given by the sublead- 
ing operators may be important and cannot be 



neglected. Using the results from steps 1 to 
4, we can evaluate the i?-parameters of all op- 
erators relevant to phenomenology. For exam- 
ple, for the renormalization scale independent B— 
parameter of Oll we found Bb^ = 1.14±0.16 and 
BbJBb^ = 1.01 ±0.01, 0. The calculation of 
the other i?-parameters is in progress and will be 
published elsewhere [||. 
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